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ABSTRACT 

We present a combined X-ray and optical analysis of the cold front cluster Abell 1201 using archival 
Chandra data and multi-object spectroscopy taken with the 3.9m Anglo Australian and 6.5m Multiple 
Mirror Telescopes. This paper represents the first in a series presenting a study of a sample of cold front 
clusters selected from the Chandra archives with the aim of relating cold fronts to merger activity, 
understanding the dynamics of mergers and their effect on the cluster constituents. The Chandra 
X-ray imagery of Abell 1201 reveals two conspicuous surface brightness discontinuities, which are 
shown to be cold fronts, and a remnant core structure. Temperature maps reveal a complex multi- 
phase temperature structure with regions of hot gas interspersed with fingers of cold gas. Our optical 
analysis is based on a sample of 321 confirmed members, whose mean redshift is z = 0.1673±0.0002 and 
velocity dispersion is 778±36 km s~^. We search for dynamical substructure and find clear evidence for 
multiple localized velocity substructures coincident with over-densities in the galaxy surface density. 
Most notably, we find structure coincident with the remnant X-ray core. Despite the clear evidence 
for dynamical activity, we find the peculiar velocity distribution does not deviate significantly from 
Gaussian. We apply two-body dynamical analyses in order to assess which of the substructures are 
bound, and thus dynamically important in terms of the cluster merger history. We propose that the 
cold fronts in Abell 1201 are a consequence of its merger with a smaller subunit, which has induced 
gas motions that gave rise to 'sloshing' cold fronts. Abell 1201 illustrates the value of combining 
multi-wavelength data and multiple substructure detection techniques when attempting to ascertain 
the dynamical state of a cluster. 

Subject headings: galaxies: clusters: individual (Abell 1201) — X-rays: galaxies: clusters 
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1. INTRODUCTION 

Within the current cosmological paradigm, large scale 
structure in the Universe is expected to form in a hierar- 
chical manner. This "bottom up" formation scenario cul- 
minates with the formation of clusters of galaxies, which 
are the largest and most massive virialized objects in the 
Universe. At the current epoch, a high fraction of clus- 
ters are still growing through the infall of matter, much 
of which is funneled through the surrounding spider-web 
like filamentary structures. The most extreme growth 
event occurs when two clusters of roughly equal mass 
merge, the merging process being one of the most ener- 
getic in the Universe, releasing around 10^'^ ergs of grav- 
itational binding energy. Some 10% of this is dissipated 
via mechanisms such as shock and adiabatic heating of 
the intracluster medium (ICM) , acceleration of relativis- 
tic particles, and generation of peculiar velocities and 
turbulence i n the ICM (fo r a review of these physical 
processes see lSarazinl[200l . 

In this context, the new generation of X-ray obser- 
vatories (XMM-Newton and Chandra) have provided a 
wealth of new information on the ICM and hence new 
insights into this cluster merger phenomenon. One of 
the first interesting discoveries - attributed to the ex- 
cellent spatial resolution and sensitivity of Chandra - 
was the observation of extended 'edge' features in the 
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X-ray surface brightnes s maps obtained for Abell 2142 
(iMarkevitch et aLlboOOl ) and Abell 3667 (jVikhlinin et all 
l2001h . The sharp discontinuitv in surface brightness seen 
in Abell 3667 had been observed previo usly with ROSAT, 
and was interpreted as a shock front (jMarkevitch et ahl 
Il999f ). However, the Chandra observations showed that 
in both Abell 2142 and Abell 3667, the temperature of 
the gas on the brighter (denser) side of the discontinuity 
was colder than the temperature of the less dense down- 
stream gas - opposite to expectations for a shock. Fur- 
thermore, the temperature and density profiles across the 
discontinuities showed, in combination, that the pressure 
is roughly continuous, again inconsistent with a shock 
front. Thus, the observed edges appear to be contact 
discontinuities between cool, dense, low entropy gas and 
hotter, diffuse ambient ICM, leading them to be dubbed 
"cold fronts" . Initial explanations attributed the low en- 
tropy g as to the remnant coolin g core of a merging sub- 
cluster (jMarkevitch et al.ll2000l ). 

As the number of clusters observed to have cold front 
features has increased, it has become clear that the rem- 
nant core scenario is not the correct interpretation in 
all cases, since a number of cold fronts are observed 
to exist in clusters with an other wi se relaxed X-ray 
morphology (Mazzotta et al.' '20013,"^; IMarkevitch et all 
120011: IMarkev itch & Vikhhnin 2007,). One such ex- 
ample is Abell 1 795, the observations of which led 
IMarkevitch et al.l ()2001l ) to propose an alternative sce- 
nario whereby the 'sloshing' of the cooling core within 
the stationary gravitational potential well causes cool 
central gas to be displaced and a cold front formed where 
it comes into contact with higher entropy gas at larger 
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radii. 

Hydrodynamic simulations have an integral role in aid- 
ing the interpretation of these cluster X-ray observations 
and further understandin g the under lying physics . Both 
lAscasibar fc Markev itch (2006*) and iPoole et all (120061 ) 
found a plethora of transient cold-front like phenom- 
ena associated with sub-cluster gas arising during their 
simulated mergers, including the classic ram-pressure 
strippe d remnant sub-cluster core p receded by a merger 
shock. iTittlev fc HenriksenI (|200 5^ explored the possi- 
bility of forming a cold front through gas sloshing, find- 
ing it is possible to produce similar edges to those ob- 
served through core oscillations, although in their model 
they find the gravitational potential well oscillates due 
to the motion of the da r k ma tter, inducing motion in 
the gas. iChurazov et all (|2003f ) and iFuiita e t al. (200|) 
also showed that gas oscillations can be induced when 
the cluster core is displaced from the potential well by a 
weak shock or ac oustic wave which has passed t hrough 
the cluster center. lAscasibar fc MarkevitchI (|2006f ) found 
that oscillation of the central gas and dark matter can 
easily be induced by infalling sub-clusters, and can even 
be induced due to the infall of a dark matter only sub- 
cluster. The gas core decouples from the dark matter 
when a rapid change in the direction of oscillation due 
to the core passage of the dark matter sub-cluster causes 
a change in the ram pressure felt by the gas, displacing 
it from the dark matter. The gas sloshing then occurs 
when the displaced gas falls back toward the potential 
minimum, generating edges at the turn-around point of 
each oscillation. 

Whatever the mechanism, it appears the existence of 
a cold front can generally be interpreted as strong ev- 
idence of a system which is in the process of, or has 
recently undergone a merger. What has been lacking 
to date is a systematic study of the relationship be- 
tween cold fronts and other dynamical indicators of clus- 
ter merger activity using multi- wavelength observations, 
which can be incorporated into detailed models in or- 
der to garner a complete understanding of how a ma- 
jor merger impacts the different cluster mass compo- 
nents. To address this issue, we have conducted a search 
of the Chandra archives and selected a sample of clus- 
ters exhibiting robust examples of cold fronts (the se- 
lection criteria and sample will be presented in a forth- 
coming paper; Owers et al. 2008, in prep.) for optical 
multiple-object spectra (MOS) and radio follow up obser- 
vations. This paper represents the first of a series where 
we aim to show a relationship between cold fronts and 
merger activity through detection of substructure at op- 
tical wavelengths. Subsequent papers will relate the star 
formation/radio properties of the galaxies and large scale 
diffuse r adio halo/relic emission ass ociated with cluster 
mergers (jGiovannini fc Ferettill20f)l . 

Detection of substructure within clusters using optical 
data has a long history. Initially, galaxy surface density 
contours were used to searc h for projected galax y con- 
centrations within clusters (jGeller fc Beer£iri982f ). Ad- 
vances in multi-object spectroscopy have allowed simul- 
taneous observations of tens, and now hundreds of galaxy 
spectra within a field, meaning large samples of spectro- 
scopically confirmed cluster members can now be com- 
piled. Radial velocity information aids in eliminating 
the projection effects inherent in galaxy surface density 



contour methods, and a number of methods have been 
proposed which use the combination of galaxy radial ve- 
locity and spatial information to develop a statistic for 
the robust detection of sub structu res within a cluster (eg . 
see Dressier & Shectman 1988; West fc Bothun 199QI; 
[CoUess fc Dunn i9_96,; Girardi fc Biviano ,2002). Since 
different s tatistics are sensitive to different types of sub- 
structure (jPinknev et al.|[l996l ). the use of a combination 
of a number of different statistical methods is essential 
in diagnosing the dynamical activity in a cluster. Com- 
bining these substructure detection methods with data at 
X-ray wavelengths has proven to be a powerful tool in di- 
agnosing cluster merger scenarios and is essential for dis- 
entangling the comp l ex histories of merging systems (eg. 
Boschin et al. 2004; Girardi et al. 2006; Barrena _et al.l 
2007t ICarra sco et al. 2007; Maurogordato et al. 2006 



Here, in this first paper, we present a detailed optical 
and X-ray analysis of the cluster Abell 1201. The (previ- 
ously unpublished) Chandra observations of this cluster 
reveal two sharp surface brightness discontinuities, both 
of which are cold fronts, and also an offset core of X-ray 
emission. Based on co mparisons of the simulat i ons o f 
IPoole et al.l (|2006l ) and lAscasibar fc MarkevitchI (j2006l ) 
and these observations, Abell 1201 appears to be an ex- 
cellent example of cold fronts generated by core gas mo- 
tions caused by a gravitational perturbation in the form 
of a merging subcluster. Thus, Abell 1201 provides a 
unique opportunity to test this scenario using combined 
optical and X-ray analyses. 

Abell 1201 is a richness class 2 clu ster dA bcll et aLl 
19891) at moderate redshift z=0.168 (|Strubl c fc Roo3 



19991) with an X- ray luminosity ~ 2.4 
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(jBohringer et al.l[2000) . The analysis presented is based 
on comprehensive optical spectroscopy obtained us- 
ing both AAOmega on the Anglo- Australian Telescope 
(AAT) and Hectospec on the Multiple-Mirror Telescope 
(MMT) combined with Chandra archival X-ray data. 
The structure of the paper is as follows: In ^we present 
the reduction and analysis of the Chandra data. In ^ 
we present the new MOS optical observations and anal- 
ysis. In 21 present a scenario for the merger history 
and formation of the cold fronts in Abell 1201. In fj5]we 
summarize our results. 

Throughout the paper, we assume a standard ACDM 
cosmology where i?o = TOkms""'^, Vim = 0.3 and SIa — 
0.7. For this cosmology and at the redshift of the cluster 
1" = 2.88 kpc. 

2. X-RAY OBSERVATIONS AND ANALYSIS 

Our X-ray analysis of Abell 1201 uses Chandra archival 
data that were taken on 2004 November 3 (Obsid 4216). 
The observations were made with the ACIS-S array with 
the cluster centered on the back illuminated S3 chip, for 
a total exposure time of 40 ksec. 

2.1. Data preparation 

The data were reprocessed using the CIAO software 
package (version 3.4) starting from the standard Chandra 
pipeline processed level-1 event list. Observation specific 
bad pixel files containing hot pixel and cosmic ray af- 
terglow information were produced and applied, the lat- 
est gain files and calibrations were applied and VFAINT 
mode cleaning was used for improved rejection of cosmic 
ray events. The data were then filtered to include only 
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events with ASCA grades 0, 2, 3, 4 and 6. 

Since the backside illuminated chips are prone to flare 
contamination, we filtered the data for periods of anoma- 
lously high background. We followed standard proce- 
dure and extracted lightcurves from source-free regions 
to search for flares. Cluster emission covers much of 
the S3 chip, so we extracted a lightcurve in the 2.5- 
6keV range using the SI chip which has very similar 
flare properties to the S3 chip. The observation suffered 
from a strong flare, and approximately half the expo- 
sure time was affected, leaving a useful exposure time of 
21.5 ks. During the spectral analyses presented below, 
this cleaned data is used. For the imaging analyses, in- 
cluding the fitting of surface brightness profiles, the full 
40ks uncleaned exposure is used. 

Background data were taken from the blank sky obser- 
vations appropriate for the epoch of observation^. The 
backgrounds are taken from observations with low Galac- 
tic foregrounds and soft X-ray brightness, so we check 
the soft X-ray flux in the vicinity of Abell 1201 using 
the ROSAT all sky R45^ count rates and confirm they 
are consistent with the blank sky background rates. The 
background and source data were processed using the 
same calibration files, bad pixel files and background fil- 
tering, with the backgrounds being reprojected to match 
the observations. 

We check for residual backgrounds in the cleaned data 
by extracting spectra from the SI chip and after sub- 
tracting backgrounds from the blank sky observations, 
find a small residual background. This residual was mod- 
eled by a soft thermal MEKAL component with solar 
abundance and kT=0. 27keV (similar to that found in 
iMarkevitch et alll2003l ) plus a cutoff power law, which is 
not folded through the instrument response, with pho- 
ton index of -0-15 and exponential cutoff at 5.6 keV 
(jMarkcvit ch et aDl2003l ). Including these components, 
scaled for region size, in the spectral fits performed in 
§ 12.31 makes no significant difference to the results. 

2.2. X-ray images and Surface Brightness Modeling 

The raw 0.5-7.0 keV Chandra image, uncorrected for 
exposure and showing the entire ACIS-S3 chip, is pre- 
sented in the left panel of Figure [1] The right panel 
of Figure [T] shows the corresponding exposure corrected, 
adaptively smoothed Chandra image. The cluster has a 
bright core, and appears elongated along an axis point- 
ing to the north-west. Along this axis, there are three 
distinct features apart from the cool core: There are two 
surface brightness discontinuities, one to the south-east 
~ 300 kpc from the core running from position angle (PA; 
measured from due west) 240° to 267° and one ^ 50 kpc 
north-west of the core (PA 8° to 120°) and ~ 430 kpc 
to the north-west there is a faint, diffuse excess clump of 
emission. We discuss the discontinuities further in §2.3.21 
In Figure [2] we show the SDSS r-band optical image of 
the central regions of Abell 1201 with Chandra X-ray 
contours overlaid. The X-ray peak is offset by ~ 2.8" 
from the brightest cluster galaxy (BCG) and the excess 
clump of emission approximately coincides with a second 
concentration of galaxies, which will be discussed further 
in 33] and SI 



2.2.1. Beta model 

In order to emphasize the structures seen in F i gure [p 
we use the method of [Neumann fc Bohringeil (|1997l ). 
The surface brightness distribution is fitted with an az- 
imuthally symmetric model and the residuals are in- 
spected for significant deviations. We use the Sherpa 
fitting package to fit a background subtracted double- 
Beta model of the form 

-3/3+1/2 



5(r)=^5o,. 



1 



B, 



(1) 



^ See 
5 See 



http://cxc.haxvaxd.edu/contrib/maxim/acisbg/ 1 



htt p : / / heasarc .gstc . nasa.gov/docs / tools . html | 



where r is the radius, measured from an adjustable ori- 
gin. One Beta model accounts for the excess emission 
within the central 100 kpc and the second models the 
remaining cluster emission. For fitting, the model is mul- 
tiplied by the exposure map to account for instrumental 
effects. The best fitting model parameters are presented 
in Table [11 

We obtain residual maps by subtracting the smooth 
model described above from the entire cluster image. 
This residual map is smoothed with a Gaussian kernel 
with a = 5.9" and residual significance is deterrn ined us- 
ing the method of [Neumann &: Bohringeil (|1997( ). where 
Gaussian smoothing and Poissonian statistics are used 
to determine an error map for the X-ray image, and the 
significance map is the ratio of the residual to error map. 
Figure [3] shows the residual significance map, with the 
north-west clump detected at greater than 10a signifi- 
cance while an excess plume is detected in the direction 
of the south-east discontinuity. There appears to be a sig- 
nificant excess detected in the south-west at the bottom 
of the chip. This can be attributed to excess emission 
due to the background flare, since the excess is much less 
significant when the same analysis is performed on the 
cleaned 21.5 ks exposure. Also apparent are the signifi- 
cant negative residuals surrounding the south-east plume 
and in the outer parts of the chip to the west, artifacts 
of the azimuthally symmetric model. 

2.3. Clobal spectral properties and the nature of the 
surface brightness discontinuities 

In this section we derive global properties for 
Abell 1201, and also analyze the properties of the sharp 
surface brightness discontinuities seen in Figure [1] with 
the aim of determining the nature of the fronts. Table [2] 
summarizes the various best fitting spectral parameters 
measured in the following sections. 

2.3.1. Global Temperature 

We derived a global temperature and metallicity for 
Abell 1201 within an elliptical region covering the major- 
ity of the cluster emission (see Figure S]) . The spectrum 
was extracted using the CIAO dmextract tool, and point 
sources detected with wavdetect excised. The CIAO tool 
mkwarf was used to create the auxiliary response file 
(ARF), which accounts for spatial variations in quan- 
tum efficiency (QE), effective area and the buildup of 
contaminant on the filter windows. The observation was 
performed at — 120°C, so the CIAO tool mkacisrmf was 
used to create redistribution matrix files (RMF). The 
ARF and RMF calculated within a region are weighted 
by the number of events in the 0.5-2 keV range within 
the region. 
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Fig. 1. — Left: Raw Chandra ACIS-S3 X-ray image of Abell 1201 in the 0.5-7.0keV energy range, binned to 4" pixels, including only the 
S3 chip, using the full 40ks exposure and displayed on a logarithmic scale to enhance the diffuse, low surface brightness emission. North is 
up and east is to the left. Right: Adaptively smoothed, exposure corrected close up version of the left image with point sources removed. 
The emission appears elongated on an axis from the south-east to north-west, and there is a bright central X-ray core. To the north-west 
there is a clump of excess emission and there are two surface brightness discontinuities, one to the south-east and one near the core. 



TABLE 1 

Beta model parameters derived from fitting the surface brightness distribution 
with a double beta model (see text). 



xo,yo So,l So.2 P 1-0,1 r-0.2 B 

(deg, J2000) (10"^) (10-'^) (kpc) (kpc) (lO"") 

ll''12'"54.3=, -|-13°26'06.7" 3.53t{j|^ 0.44^0.04 O.SOtJJ IJ^ 22^^ 213^^^ 6.56tg; 



Note. — Units of -So,! , 5*0,2 and B are photons / cm^ / s/ arcsec^ . Errors on xq and yo arc ~ 0.4" . 



TABLE 2 

Summary of best fitting spectral parameters 
for several regions. 



Region 


kT 


Abundance 




(keV) 


(Z) 


Global 


5.3 ±0.3 


0.34 ±0.10 


Inside south-east front 


3.6tJ;? 




Outside south-east front 






Inside north-west front 


o 9 + 1.0 
•^■^-0.7 




Outside north-west front 


''•^-l.O 





The unbinned spectra were fitted in the 0.5-9.8 keV 
range using the MEKAL (|Kaastral Il992l : 'Liedah l et"an 
[l995h and WA BS models within the XSPEC package 
( ArnaudI [1996), minimizing the Cash statistic. The 
MEKAL component models a hot, diffuse, single tem- 
perature plasma, and is multiplied by WABS, a pho- 
toelectric absorption model accounting for Galactic ab- 
sorption, where the neutral hydrogen column density 
is fixed to the Galactic value of Nh — 1-61 x 10^" 
(|Dickev fc LockmanI 1199(1 ). The best fitting values for 
the temperatures and abundance are kT = 5.3^03 and 
Z = 0.34j^Q'|^Q where the errors are 90% confidence lim- 
its. We used the LUMIN function to determine the 
integrated unabsorbed X-ray luminosity from the best 



fitting model in the energy range 0.5-7.0 keV, and ob- 
tained Lx(0.5 - 7.0 keV) = 3.7 x 10'''' ergs"! within the 
r ~ 550 kpc elliptical aperture shown in Figure SI 

2.3.2. Surface Brightness Discontinuities 

To characterize the surface brightness discontinuities 
as either shock or cold fronts, we fitted the surface bright- 
ness profiles across them with a density model, and also 
measured the temperature difference across the discon- 
tinuities. The density model consists of a broken power 
law function such that 

I "e,2(ijj , r>Rf, 

where Rf is the radius at which the discontinuity occurs, 
rie,! and ne,2 are the densities of the inner and outer 
gas at Rf, respectively, and the spheroidal radius, r, is 
defined by = tjP' + e^^C^- Here C is the coordinate 
along our line of sight and the elliptical radius, ro, is 
defined by 

2 e^[a;cos0 + ysin6']^ + [ycos0 — xsin^)]^ 

= , (3) 

where x and y are cartesian coordinates centered at the 
center of curvature of the discontinuity and determines 
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Fig. 2. — SDSS r-band image of the central region of Abell 1201. 
The red contours are the Chandra X-ray contours logarithmically 
spaced by a factor of 1.3 in the interval 1.5 X 10~** — 4.6 X 10^^ pho- 
tons cm~^ s~^ pixel"^ , the black squares identify galaxies which 
are cluster members, the green squares show galaxies allocated to 
KMMl, and the green elli pse sh ows the Icr region for KMMl from 
the KMM analysis (see § 13. 3.31 1 . The bright cluster galaxy white 
square appears to be associated with the bright X-ray core, al- 
though is offset by ~ 2.8". 




Fig. 4. — Close up of the raw Chandra image with regions show- 
ing where surface brightness profiles were fitted, and where spectra 
were extracted for temperature measurements. The large green 
ellipse shows the region where spectra were extracted for global 
temperature and abundance measurements. The surface brightness 
profiles shown in Figure [5] are taken from the solid red and white 
sectors. The boundaries indicate the fronts, and delineate regions 
where the temperature inside and outside the fronts are measured. 
The dashed red sector indicates the slightly larger opening angle 
used in measuring the temperature outside the south-east front. 




Fig. 3. — Residual significance map derived from subtracting a 
smooth double Beta model from the X-ray image and calculating 
the significance (see text). Grey scales show the negative residuals, 
whilst color scales show the positive residuals. The black contours 
range from 1 — 12cr with linear increments of Icr. 

the orientation of the front on the plane of the sky. We 
assumed the gas density distributions inside and outside 
the front have the same eUipticities and have rotational 
symmetry so that — 1/e. The model is incorporated 
into Sherpa where we fitted the surface brightness across 



the discontinuity in the sectors shown in Figure 2] using 
the full 40 ks exposure, in the energy range 0.5-7 keV, and 
also add to the model a constant background component 
which was measured in § 12.2.11 We present the surface 
brightness distribution, along with the best fitting den- 
sity models, in Figure[5l From the amplitudes of the best 
fitting surface brightness model, Ai and we derived 
a density jump of ^J~A\fA^ — Ue^ijuej. — 2.13iQ;3g ^^r 

the south-east discontinuity and ne,i/ne,2 — 1.751q;34 
for the north-west discontinuity. The confidence range 
for the density jump was computed from the extremes 
of the 90% confidence ranges for Ax and A^, assuming a 
single parameter of interest. We note that if the density 
jump is the only parameter of interest, this overestimates 
its confidence range. 

The temperature profile across the front is an essential 
diagnostic tool. Ideally, to minimize projection effects we 
would like to measure the deprojected temperature pro- 
file across the fronts, however the limited number of pho- 
tons available do not allow this. Instead, we extracted 
spectra from two regions for each front, one on the bright 
side of the front (inside) and one on the faint side (out- 
side). The regions from which we extracted spectra for 
both the south-east and north-west fronts are shown in 
Figure m For the north-west front, we extracted spectra 
from one sector centered at RA=168.2268, DEC=13.4344 
with PA from 7° - 122° and radii - 18.4" and one with 
the same center and PA, but with radii 18.4 — 45.5", cor- 
responding to the sector in which we fit the density model 
above. For the south-east front, we extracted spectra 
in one sector centered at RA=168.2197, DEC=13.45094 
with PA from 240° - 267° and radii 127.9 - 175.1", and 
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Fig. 5. — Left: Surface brightness profile across the north-west discontinuity measured in a sector from PA 7° — 122° centered at the center 
of curvature of the front (RA=168.2268, DEC=13.4344).Ri5/ii; Surface brightness profile across the south-east discontinuity measured in 
a sector from PA 240° — 267° centered at the center of curvature of the front (RA=168.2197, DEC=13.4509). The red curves show the 
surface brightness profile for the best fitting density model. 

the temperature on the dense side of the front) and com- 
paring with the values observed. For shock fronts, we 
would expect to observe pre-shock gas temperatures of 
1.9^0^4 keV and 2.11q g keV for the south-east and north- 
west fronts, respectively, significantly different from our 
measured temperatures. With continuous pressure and a 
temperature increase, these discontinuities have the hall- 
marks of cold fronts. 



one with the same center, but PA from 234° — 273° and 
radii 175.1 — 270.6". Note that the region from which we 
extracted a spectrum for outside the south-east front has 
a slightly larger opening angle compared to the region in 
which the surface brightness profile was fitted. The sec- 
tor is wider for better statistical accuracy. Spectra and 
responses for these regions were extracted as described 
in g2Xll 

The spectra outside and inside the fronts were fitted 
successively, with the regions outside the front fitted us- 
ing an absorbed MEKAL model with the column density 
and metallicity fixed to the values derived in ^12.3.11 The 
regions inside the fronts were then fitted with a model 
containing two absorbed MEKAL components, the pri- 
mary MEKAL component models the gas within the 
front and the secondary MEKAL component accounts 
for gas lying in projection along the line of sight, which 
is assumed to have the same thermal properties as the 
gas lying in the region outside the front. The column 
density and abundance were fixed to the values derived 
in ^ 32.3.11 for both components and the secondary com- 
ponent had its temperature fixed to the best fit value 
obtained from the outer region spectra. The normaliza- 
tion of the secondary component was fixed to the value 
obtained from fitting the spectrum outside the front, cor- 
rected by a factor that accounts for the different emission 
measures expected from the different volumes probed, 
which is calculated by integration of the density model. 
For the south-east front, we measured a temperature of 
5.7l^'7keV outside the front and 3.6^QjkeV inside the 
front, and for the north-west front we measured a tem- 
perature of 5.2t_i^keV outside the front and 3.2^0.7 keV 
inside the front. 

The pressure should be continuous at a cold front. The 
pressure jump across a front can be measured by tak- 
ing the ratio {rie^ikTi) / {ne,2kT2). We measured pressure 
jumps of 1.4]';|!) g and I.IIq g for the south-east and north- 
west fronts, respectively, consistent with the pressure be- 
ing continuous across both fronts. A shock interpreta- 
tion for the fronts can be also ruled out by applying the 
Rankine-Hugoniot shock jump conditions for the mea- 
sured density jumps and post-shock temperatures (i.e. 



2.3.3. Temperature Map 

Temperature maps provide extremely useful tools for 
searching for multi-phase temperature structure due 
to the effects of an ongoing merger. To search for 
evidence of merger induced temperature structure in 
the ICM of Abell 1201, we generated a temperature 
map using the br o ad en ergy band method described in 
iMarkevitch et all (|2000r ). Briefly, we produced source 
and background images, binned in 7.8" pixels, in 
the energy bands 0.5-1. 0-2. 0-5. 0-lOkeV, excluding point 
sources. Exposure maps that corrected for mirror vi- 
gnetting relative to the on-axis position, QEU (includ- 
ing low energy contamination) and exposure time were 
produced for each energy band. The background im- 
ages, which were taken from the blank sky observations 
described above, were normalized by the ratio of the 
source to background 9-12 keV counts, subtracted from 
the source images and the resulting image divided by its 
corresponding exposure map. Each image was smoothed 
using the same variable width Gaussian, where a var- 
ied from 8.6-39.4", with it being smallest in the bright- 
est regions and chosen such that statistically significant 
temperatures could be measured whilst maximum spatial 
information was retained. The noise in each pixel was de- 
termined from the raw, uncorrected images and weighted 
accordingly to allow for the effects of the smoothing. 

Each pixel was fitted with an absorbed single temper- 
ature MEKAL model with the absorption column set at 
the Galactic value. The metal abundance was set to 
the best fitting average cluster value derived in i j2.3.1l 
The model was multiplied by an on-axis ARE to cor- 
rect for the energy dependent on-axis mirror effective 
area, including the chip quantum efficiency (QE). A flux- 
weighted spectral response matrix was generated from a 
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Fig. 6. — Temperature map with 0.5-7.0 keV brightness contours 
overlaid (same as Figure [2]|. The 68% temperature uncertainties 
range from ±0.4 keV in the cooler (kT < 5keV) regions to ±1 — 2 
in the hotter (kT > 5 keV) regions. The color differences represent 
regions which are at significantly different temperatures. The color 
bar shows the temperature scale in keV. 

large cluster region, and was binned to match the chosen 
energy bands. The temperature map is presented in Fig- 
ure [6l where we have excluded pixels where the la error 
is greater than 30% of the best fit temperature value. 

The cold fronts described above are clearly visible in 
Figure [5] and the temperature map values are consistent 
with the temperatures measured spectroscopically in the 
above sections. Interestingly, there is a finger of cool gas 
extending from the south-east cold front to the core of 
the cluster. This finger is coincident with the residual 
observed after subtraction of the Beta model in §2.2.11 
Also of interest is the significantly hotter ~ 7.5 keV re- 
gion lying between the core and the excess subclump. 
The temperature in the region containing the north-west 
excess appears to be consistent with the global cluster 
temperature. Since the surface brightness here is low, 
the photon numbers are also low, meaning the region 
must be heavily smoothed to obtain statistically signif- 
icant temperatures (u ~ 20" here). Thus the temper- 
ature measurement is averaged over a large region and 
is heavily contaminated by ambient gas projected onto 
this region, thus the temperature of the excess remains 
uncertain. 

3. OPTICAL SPECTROSCOPY - REDUCTION AND 
ANALYSIS 

In this section we present our optical analysis which in- 
cludes selection of photometric samples for spectroscopic 
follow up, spectroscopic observations and data analysis, 
determination of cluster membership and substructure 
detection. The purpose of this optical analysis is to 
search for substructures and use them in correlation with 
the X-ray observations to develop a scenario for the for- 
mation of the cold fronts observed above, and also for 
the merger history of Abell 1201. 



The multi-object spectroscopic (MOS) data presented 
in this paper were taken from two sets of observations, 
one set at the 3.9m Anglo-Australian Telescope (AAT) 
on 2006 April 2-5, and the second at the 6.5m Multi- 
ple Mirror Telescope (MMT) in queue schedule mode in 
the months of February and April 2007. The observa- 
tion details are hsted in Table [3] where we list the dates, 
magnitude limits, frame exposures and the seeing. 

3.1.1. Parent Photometric Catalog 

The initial target catalog was taken from the Sloan 
Digital Sky Survey (SDSS) sky server^, with all objects 
within an ISarcmin radius of the center of Abell 1201 
(RA=llhl2'"54.5", DEC=-H13°26'09.0") included. This 
catalog was then filtered to include only those objects 
classified by the SDSS pipchne as galaxies (SDSS class 3). 
The SDSS r and g magnitudes for the remaining objects 
were then converted to Jo hnson B and R magnitudes 
via equations A5 and A7 of ICross et al.l ([2004) and only 
those galaxies with R < 21.5 were retained in the final 
target Hst. 

The target galaxies were then ranked based on their 
cluster-centric radius and their position on the B-R vs 
R color-magnitude (CM) diagram, in preparation for the 
fiber allocation procedure. Galaxies on or blueward of 
the red sequence are more likely to be cluster members 
and so were ranked higher than those redward of the red 
sequence. Those galaxies either on or blueward of the red 
sequence were then ranked by cluster-centric distance, 
with those closest to the center ranked highest. The re- 
maining galaxies lying redward of the red sequence were 
also ranked by cluster-centric distance. These rankings 
were input into the software used to configure the fiber 
allocations for the respective observations, to ensure that 
higher weightings were placed on galaxies which were 
more likely to be cluster members (i.e. those on and 
blueward of the red sequence and closer to the center). 
We discuss the configuration process further below. 

We show the CM diagram in Figure [7l where we have 
overplotted the spectroscopically confirmed cluster mem- 
bers (green stars), foreground galaxies (blue squares) and 
background galaxies (red circles). Also plotted is the line 
B ~ R — 3.6 — 0.0645i? which was used to distinguish 
those galaxies which we define to be lying on or below 
the red sequence from those redward of the red sequence. 
The slope of the line is calculated from the best fit to the 
CM s lope versus redshift diagram of iLopez-Cruz et all 
(l 2004l ) and the constant is estimated by eye such that 
the line lies just above the red sequence. 

3.1.2. AAT AAOmega Observations 

The AAT observatio ns were taken using the AAOmega 
fiber-fed spectro graph (' Saunders et al.|[2"004HSmith et al.l 
l200llSharD eT al. 2006), which is a bench mounted dual- 
beam spectrograph, which is fed by 400 fibers robotically- 
placed within the two degree field at the telescope's prime 
focus. A total of 392 fibers, each 2" in diameter, are 
available for the simultaneous observation of scientific 
targets (with the remaining 8 fibers being used for ac- 
quisition and guiding). All our observations were taken 
using the medium resolution (R ~ 1300) 580V (blue arm) 
and 385R (red arm) gratings, which in combination with 
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TABLE 3 



Summary of the observations. 


_LelescopG/ instrument 




iVlclgnitUQG 


r rjimcs 


■ 

oeeing 


AAT/AAOmega 


2006 Apr 2 


19.5 <R< 20.5 


3 X 1500S+ 


1.5-2.0" 








3 X 750s 






2006 Apr 4 


19.5 <R< 20.5 


3 X 1800s 


2.0" 




2006 Apr 5 


R < 19.5 


3 X 1200s 


2.0" 




2006 Apr 5 


R < 19.5 


3 X 1200s 


2.0" 




2006 Apr 5 


R < 19.5 


3 X 1200s 


2.5" 


MMT/Hectospec 


2007 Fob 21 


20.5 < R< 21.5 


4 X 1800s 


0.6" 




2007 Feb 22 


20 < < 20.5 


4 X 1200s 


1.2" 




2007 Apr 18 


R < 20.5 


3 X 1200s 


1.2" 




2007 Apr 19 


R < 20.5 


4 X 1200s 


0.8" 



R magnitude 

Fig. 7. — Color-magnitude diagram for all galaxies within our 
observed field (black dots). The green stars represent cluster mem- 
bers, the blue squares foreground galaxies, and the red circles back- 
ground galaxies. The pink line is the line used to delineate galaxies 
on and blueward of the red sequence from those redward of the red 
sequence for the purpose of ranking galaxies during fiber configu- 
ration. 



the 15 /im pixel 2kx4k E2V CCD detectors produce a 
spectral resolution of 3.6A and 5.5A in the blue and 
red arms, respectively, and provide an overall wavelength 
coverage ranging from 3700 to 88OOA. 

Due to the clustered nature of the galaxies, fiber po- 
sitioning constraints (minimum separation of ~ 30" due 
to the physical size of the fiber buttons) and the fact 
our input catalog contained ~ 1400 galaxies, multiple 
plate configurations were required to obtain high over- 
all spectral completeness and adequate coverage in the 
central regions of the cluster. The fibers were allocated 
using the AAOmega CONFIGURE software'^. A total 
of 5 configurations were required for Abell 1201: three 
for galaxies with R < 19.5, and two for galaxies in the 
range 19.5 < R < 20.5 plus those that remained unallo- 
cated in the three 'bright' configurations. Due to packing 
constraints, only 110-210 of the available fibers were al- 
located to science objects. In addition, 30-40 fibers were 
allocated to blank sky regions for sky subtraction. All 
fibers known to suffer from interference fringing effects 
were left unallocated. The details for each configura- 
tion observation are listed in Table O The observations 
were taken in rather mediocre conditions, with the see- 
ing (FWHM) ranging from 1.5 to 2.5arcsec, and all ex- 



See: |http://www. aao.gov.au/AAO/2df71 



posures being taken through thin cloud cover, with some 
being interrupted by thicker clouds. In commencing the 
observations for each new configuration, a 4s dome flat 
and 60 s FeAr arc lamp exposure were taken for flat field- 
ing and wavelength calibration of the data, respectively. 
The data were reduced with the AAO 2dFDR pipeline 
software with a patch included to fix fiber to slit position 
mapping errors present in AAOmega data taken before 
August 2006. 

Redshift identification and measurement for each spec- 
trum was carried out using the RUNZ code written by 
Will Suth erland for the 2dF Galaxy Redshift Survey 
(2dFGRS; Colless et al.ll200lD. This progra m utilizes the 
cross-correlation method of' Tonrv fc DaAdi (fT979l) . based 
on a library of galaxy template spectra that are rep- 
resentative of all the different observed spectral types. 
Each spectrum was inspected visually and given a red- 
shift quality classification, Q. Here we used the same 
scheme as adopted for the 2dFGRS, with each spectrum 
being assigned a Q value on a six-point integer scale, with 
Q=l indicating that no redshift could be estimated, Q=2 
a possible but unreliable redshift, Q=3 a probable red- 
shift (with ~ 90% confidence), Q=4 a reliable redshift, 
Q=5 a reliable redshift with high-quality spectrum, and 
Q=6 indicating a star or non-extragalactic object. We 
obtained spectra for 917 galaxies during the run, which 
yielded reliable (Q=3, 4, or 5) redshift measurements for 
580 single galaxies. 



3.1.3. MMT Hectospec Observations 

The MMT observations were taken using the Hec - 
tospec multi-object spectrograph (|Fabricant et al.|[2005f) . 
This is also a bench-mounted spectrograph, which is fed 
by 300 1.5" diameter fibers that cover a 1 degree field 
of view. The observations were taken using the medium 
resolution (R=1000-2000) 270 groove mm~^ grating and 
the data were captured on a single array of two E2V 
CCDs with 13.5/iTO pixels, resulting in spectra with a 
resolution of 6.2A and covering the wavelength range 
3500- lOOOOA. 

For the MMT observations, we included any galaxy not 
observed during the AAT run, as well as galaxies that 
were not assigned reliable redshifts. A small number of 
galaxies which were assigned a redshift quality of Q=3 
(from their AAT spectra) and were in the redshift range 
0.15-0.18 (i.e. close to the cluster in redshift space), were 
also included in order to check on our redshift accuracy. 
Fiber configurations were generated using the XFITFIBS 
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software , with a total of four being required for galaxies 
with R < 20.5, and one being required for galaxies with 
20.5 < R < 21.5. Typically 150 fibers were allocated 
to galaxies, with around 100 fibers allocated to blank 
sky areas. The details of each observed configuration 
are listed in Tabled The seeing during the observations 
ranged from 0.61-1.2", and the observations performed in 
the worst seeing conditions were somewhat affected with 
a number of galaxies requiring re-observation since the 
spectra were inadequate for determining redshifts. The 
data were reduced at the Telescope Data Center (TDC)^ 
using the TDC pipeline. The data were also redshiftcd 
at the T DC using the IRA F cross-correlation XCSAO 
software ([Kurtz et al.lll99^ and spectra were assigned 
a redshift quality of "Q" for reliable, "?" for question- 
able and "X" for bad redshift measurements. We use 
in our analysis only those galaxies with reliable redshift 
measurements (quality=Q). A subsample of the data was 
visually inspected and the redshifts and quality assign- 
ments were found to be robust. We obtained spectra 
for 742 galaxies during the run, which yielded reliable 
redshift measurements for 534 single galaxies. 

3.1.4. Redshift completeness and measurement errors 

A number of galaxies were observed multiple times in 
the two observing runs. There were also 67 galaxies 
with independently measured redshifts within 20 arcmin 
of the cluster center in the NASA Extragalactic Database 
(NED)^°. To check the accuracy of our measured red- 
shifts, we compared the redshifts of all multiply-observed 
galaxies by determining the mean of the differences as 
well as their standard deviation. The MMT observations 
had a mean intrinsic uncertainty of 42 ± 1 km s^ ^ for 534 
redshifts. There were 19 repeat observations made with 
the MMT, and the mean difference was 1 ± 25kms~^ 
after removing a spurious result where the difference was 
cz = 44894 km s^^. From these repeat observations, we 
calculate an RMS of 105kms~^ implying an uncertainty 
of 74kms~^ for the MMT observations, higher than the 
mean of the individual redshift measurements. An ex- 
ternal test of these redshift uncertainty measurements 
comes from the 12 redshifts in common between the 
NED and MMT catalogs which have a mean difference 
of -28 ± 20kms~^ with an RMS scatter of 74kms~^ 
The NED re dshifts come prima rily from the SDSS and 
the catalog of lMiller et al.l ([2002f ) which have redshift un- 
certainties ~ 30kms~^. Taking this into account gives 
an external uncertainty measurement of 43kms~^, con- 
sistent with the mean of the individual uncertainty mea- 
surements. The high uncertainty value given by the re- 
peat MMT observations can be understood by consider- 
ing that these objects were re-observed due to the low 
quality of their initial redshift measurements. 

There were no repeat observations taken at the AAT, 
which have a mean internal uncertainty for the individ- 
ual redshift measurements of 109 ± 6kms~^ from 580 
redshifts. An external check of this value comes from 
31 galaxies re-observed with the MMT which had both 
Q=3 redshifts from the AAT data set and quality=Q 
from the MMT. The mean difference in the measured 



® see: http:/ /cfa- www. harvard. edu/mmti/hectospec. html 
^ see: http: / /tdc- www. harvard.edu/ 1 
see: ht tp://nedwww.ipac.caltech.edu/ 1 
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AAT redshift {Orange, Light Biue) or MMT redshift (Pink, Darl^ biue) 



Fig. 8. — Redshift comparison for Abell 1201 using our MMT and 
AAT redshifts, and also a sample of NED redshifts. The compari- 
son of different measurements are plotted with the following color 
codes: MMT/MMT = pink, AAT/MMT = light blue, AAT/NED 
= orange, MMT/NED = dark blue. The black line represents the 
one-to-one relationship; the general scatter of the data about this 
line indicates the absence of any gross systematic errors in our mea- 
surements, and indicates our AAT and MMT measurements have 
a precision of 159 km s^^ and 74kms~^, respectively. We assume 
an uncertainty of 30kms~^ for the NED measurements. 

redshifts of these was 50 ± 44kms~^. Two galaxies 
with obviously spurious redshift differences of A(cz) = 
53417kms~^ and A(cz) = 119800 km s~^ were excluded 
from this calculation. The RMS scatter in the differ- 
ence was 248 km s^^, and taking into account the uncer- 
tainty derived above for the MMT redshifts (74kms^^), 
we calculate an uncertainty of 159kms~^. This value 
is higher than the mean of the individual redshift un- 
certainties, however again it is noted that these objects 
were re-observed due to their poor initial redshift mea- 
surements (all had Q=3), thus the measured uncertainty 
is expected to be higher. This value can be checked us- 
ing the 45 measurements in common with the NED cat- 
alog, where the mean difference was 33 ± 29kms^^, and 
an RMS scatter in the differences of 191 km s^^, after 
two spurious results where the redshift differences were 
A{cz) = 3480 km s"i and A{cz) = 10472 km s"^ were ex- 
cluded. Based on these redshift differences, we measure 
an uncertainty of 132kms~^, again after accounting for 
the ~ 30kms~-'^ uncertainty associated with the NED 
redshifts. 

Figure |8| shows the comparison of redshift measure- 
ments using different instruments in the redshift range of 
the cluster. The figure, in combination with the redshift 
differences above, show that there are no gross systematic 
errors in our measurements and the scatter about a 1-to- 
1 relationship is well encompassed by the uncertainties 
of 159kms~^and 74kms~ifor AAT and MMT, respec- 
tively, which we take from the more conservative higher 
values derived above from the double observations, and 
these values indicate the precision of our measurements. 

For galaxies which had multiple redshift measure- 
ments, the following approach was taken to determining 
what their final adopted redshift would be: If a galaxy 
was observed with both the MMT and AAT, we used the 
MMT redshift since the MMT redshifts were in general 
more precise with lower errors (as expected since the re- 
observed galaxies were ones with AAT quality Q=3). If a 
galaxy was observed more than once with the MMT, we 
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took the redshift based on the highest cross-correlation 
coefficient. If a galaxy was observed more than once with 
the AAT, we took the redshift with the highest Q value. 
Note that in no cases did we ever use the NED redshift 
if that galaxy had been observed on cither the AAT or 
MMT. When added to our single redshift measurements, 
this yielded a total of 560, 534 and 10 redshifts mea- 
sured with the AAT, MMT and sourced from NED, re- 
spectively, giving an overall total of 1104 robust redshifts 
acquired within the cluster field. 

We determined the spectroscopic completeness of the 
sample by measuring the ratio of galaxies in the parent 
photometric catalog which have reliable redshifts (as de- 
fined previously) to those which have no reliable redshift 
measurement. The spectroscopic completeness as a func- 
tion of cluster-centric radius for the R magnitude inter- 
vals 0-18, 18-19.5, 19.5-20.5, and 20.5-21.5 is plotted in 
Figure [9] where we also plot the spectroscopic complete- 
ness within 3.5 Mpc as a function of R magnitude. We 
achieve > 80% spectroscopic completeness at all radii, 
apart from within the 3-3.5 Mpc range, for magnitudes 
brighter than 19.5, whilst we obtain ~ 60% spectroscopic 
completeness at all radii for magnitudes between 19.5 and 
20.5 and ~ 20% spectroscopic completeness for magni- 
tudes between 20.5 and 21.5 within a radius of ~ 2.5 Mpc. 
Thus our catalog is very well sampled at all radii for 
R magnitudes brighter than 20.5, which corresponds to 
~ 2.3 mags down the cluster lumin osity function (assum- 
ing = -21.3 (|Yagi et al.ll2002f ). and neglecting galac- 
tic extinction and K-correction terms). 

3.2. Cluster member selection 

Identification of cluster members from our spectro- 
scopic redshifts was achieved through the elimination 
of foreground and background galaxies along the line 
of sight to Abell 1201. An initial rejection was per- 
formed using th e "velo city gap" method outlined by 
iDe Propris et al.) ()2002l ). where the galaxies are sorted 
in redshift space and the velocity (cz) gap between each 
one determined. Here the velocity gap for the nth galaxy 
is Avn = czn+i — czn- Clustcrs appear as well populated 
peaks in redshift space which are separated by velocity 
gaps of greater than 1000 km s~^ from the nearest fore- 
ground and background galaxies, and Figure [10] shows 
that Abell 1201 is readily identified on this basis. 

Due to the filamentary structures which surround clus- 
ters, it is necessary to further refine this rejection pro- 
cess. For this refinement, we used a slightly different ver- 
sio n of the "shi f ting g apper" method first implemented 
bv lFadda et al.l (|1996D where both cluster-centric radius 
and peculiar velocity information are used and the gap 
method outlined above is applied as a function of radius. 
The galaxies were binned radially such that each bin con- 
tained 35 galaxies, and were then sorted by Wpec, with the 
velocity gaps determined as before (but in peculiar ve- 
locity, not cz). Peculiar velocities were determined with 
respect to the central cluster redshift using the following 
procedure: We estimated the m ean cluster redshif t using 
the biweight location estimator (jBeers et al.lllQQOl ) which 
we assume represents the cosmological redshift of the 
cluster. The peculiar redshift, Zp^c, was derived under 
the assumption that the observed redshift of the galaxy, 
Zgai, is comprised of only two components, the cosmologi- 
cal component, Zcos, and the component due to the pecu- 



liar motion within the cluster, Zp^c- Hence, the peculiar 
redshift is Zp^c = [zgai — 2cos)/(l + Zcos) and the peculiar 
velocity was derived using the standard special relativis- 

tic formula Vpec = c((l -|- Zpec)^ - !)/((! + Zpec)'^ + 1), 

where c is the speed of light. 

We used the "f pseudosigma" (jBeers et al.lll990( ). de- 
rived from the first and third quartiles of the peculiar 
velocity distribution, as the fixed gap to separate the 
cluster from interlopers. The f pseudo-sigma is an esti- 
mator of the scale of the velocity distribution which is 
robust to the presence of velocity interlopers in the tails 
of the cluster distribution. The above procedure was 
iterated until the number of members was stable and 
the results are shown in Figure [TT] where it can be seen 
that the cluster is clearly separated from the filamentary 
structure surrounding it, and the interlopers are cleanly 
rejected. The advantage of this method is that it does 
not assume a particular mass model or velocity distribu- 
tion which both rely on the cluster being relaxed. 

The final cluster sample contains 321 members out to 
a cluster-centric radius of 3.5 Mpc. The final value for 
the biweight location estimator of the cluster redshift is 
Zcos — 0.1673ib0.0002 . We use the biweight scale estima- 
tor (jBeers et al.lll990D to estimate a velocity dispersion of 
778 ±36 km s~^. The errors for the redshift and velocity 
dispersion are la and are estimated using the jackknife 
resampling technique. 

3.3. Substructure Detection 

Now that we have separated our cluster members from 
foreground and background interlopers, we can apply 
substructure detection tests to the cluster member sam- 
ple. Substructure in a cluster can present itself in a num- 
ber of different ways, and generally no single statistical 
test is capable of revealing all such manifestations. Hence 
it is essential that the full range of statistical tests are ap- 
plied if the search for substructure is to be an exhaustive 
one. In general, statistics using the maximum amount of 
information (eg., the radial velocity and spatial dimen- 
sions) are the most effective at revealing substructure, 
however these can fail if, for example, roughly equal mass 
clusters are merging along the line of si ght such that 
the cores are spatially coi ncident (see eg. iGirardi et al.l 
120061: iPinknev et "all 119961 ). The use of velocity infor- 
mation by itself can be successful, although there are 
documented cases where the velocity distribution mim- 
ics that of a relaxed cluster where other methods have 
clearl y shown it to be distu rbed (jMaurogordato etldl 
120001: 1 Johnston-Hollitt efall |2008[ ) . The use of 2-D spa- 
tial information alone suffers from foreground and back- 
ground contamination. In this section, we apply a num- 
ber of statistical methods to the detection of substructure 
in Abell 1201. 

3.3.1. Deviations from Gaussianity 

The velocity distribution of a relaxed cluster is well 
approximated by a single Gaussian, with velocity dis- 
persion related to the cluster mass via the virial the- 
orem. Departures from Gaussianity can be attributed 
to a dynamically active cluster where substructures may 
cause symmetric and asymmetric distortions in the veloc- 
ity distribution. As a first test of Gaussianity, we used 
the standard Kolmogorov-Smirnov (K-S) test. At the 
90% confidence level, the observed velocity distribution 
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Fig. 9. — Spectroscopic completeness as a function of cluster-centric radius {left panel) for the ij-band magnitude intervals 0-18 {solid 
line), 18-19.5 (dashed line), 19.5-20.5 (dot-dashed line), and 20.5-21.5 (dotted line). This shows we have both good completeness and radial 
coverage for magnitudes brighter than _R = 20.5. The spectroscopic completeness as a function of magnitude is shown in the right panel. 
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Fig. 10. — A histogram of all the reliable redshifts measured 
within the Abell 1201 field. The Abell 1201 cluster clearly stands 
out against the foreground and background galaxies. The inset 
panel shows the distribution of peculiar velocities for the cluster 
members, with a Gaussian with mean zero and standard deviation 
of 778 km s~^ overplotted. 

is consistent with being drawn from a Gaussian distribu- 
tion with = and a^pec — 778 kms^^. On the basis 
of this test, therefore, it appears unhkely AbeU 1201's 
velocity distribution is significantly non-Gaussian. 

The disadvantage of using the K-S test is that it is 
most sensitive to the behavior of the distribution near its 
median, but is relatively insensitive to differences in the 
tails of the distribution. Also, the K-S test does not give 
quantitative information about the way in which two dis- 
tributions differ. Quantifying these deviations is impera- 
tive in determining the state of the cluster, and a number 
of different methods can be use d to quant ify asymmetric 
and symmetric distortions (Pin knev et a l. 1996) . Here 
we cho ose to use the method outlined in Zabl udoff et al.l 
(|1993f ). where the velocity distribution, L, is approxi- 
mated by a series of Gauss-Hermite functions 

4 

L = J2h,H,{x), (4) 
where hj are the Gauss-Hermite moments defined by 



Fig. 11. — A plot of deviations in peculiar velocities as a func- 
tion of cluster-centric radius, illustrating the efficacy of our refined 
rejection method whereby foreground and background interlopers 
are identified and eliminated using a shifting gapper technique (see 
text). The black crosses represent galaxies allocated as cluster 
members while the orange asterisks are rejected foreground and 
background galaxies lying close to the cluster in redshift space. 

IZabludoff etaH (Il993l) as 

i—l 

where 

^ V_^ec^^ ^^^^^ ^ t__n,{xl (6) 

and Ti.i(x) are the Hermit e polynomials given by 
Ivan der Marel fc Franxl (|1993D . In principle, the veloc- 
ity distribution can be described by equations [4] and 
[5] and a wide variety of S and V values, where there 
is significant degeneracy in the choices of S and V 
and the j = 1,2 terms. Ho wever, as pointed out by 
Ivan der Marel fc Franxl ()1993f l , the most efficient method 
is to set S and V so that the zeroth-ordcr term, Hq, de- 
scribes the best fit Gaussian to the data. This is achieved 
by choosing S and V such that hi = h2 = 0. We iterated 
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with different values of S and V until these criteria were 
met. The terms /13 and /14 describe, respectively, the 
asymmetric and symmetric deviations from Gaussianity, 
much like the higher order Gaussian skewness and kurto- 
sis terms, but have the advantage of being less sensitive 
to outliers in the tails of the distribution. Applying this 
method to the velocity distribution for Abell 1201, we 
obtain values of values of /13 = —0.016 and — 0.023. 
We used 10,000 Monte Carlo realizations of Gaussian dis- 
tributions with N=321, a = S and ^ = V io determine 
the probability of detecting our /13 and /14 terms. Values 
of \h^\ > 0.016 occur in 70% of the realizations, while val- 
ues of I/14I > 0.023 occur in 50% of the realizations. We 
conclude there is no significant evidence for asymmetric 
or symmetric deviations from Gaussianity. 

3.3.2. Velocity Distribution Profiles 

Figure [12] shows integral and differential projected ra- 
dial profiles of the biweight location, and scale, 
cr{R), estimators. The differential profiles are binned us- 
ing the same radial intervals as used for the shift gap- 
per in §3.2^ such that each bin contains at least 35 data 
points. To measure the integrated profiles, we sorted 
the data by cluster-centric radius and for each galaxy in 
the sample, starting from the tenth, we measured the bi- 
weight location and scale (where the number of galaxies 
was less than 15, we used the median and gapper esti- 
mates of location and scale) using only data within the 
radius of interest to each galaxy. 

The integrated and differential /i(-R) are constant with 
radius within the la error bars. This is expected for an 
isotropic distribution of radial velocities. The differential 
(j{R) profile is flat o ut to r = 1 Mpc, and smooth ly de- 
clines at larger radii. Iden Hartog fc Katgertl (|1996f l clas- 
sified (j{R) profiles as flat, inverted or peaked, depending 
on the shape of the profile within ^ 1 Mpc, and they note 
an interesting subset of clusters with flat a{R) profiles 
which have the highest cr(i?) value measured between 
0.5 and 1.0 Mpc, all of which show signs of recent merger 
activity. While this is by no means a definitive sign of 
merger activity, it is interesting to note that Abell 1201 
also exhibits its highest value of cr(i?) at ~ 1 Mpc. The 
smooth decline outside of IMpc is expected since the 
escape velocity of a cluster declines with radius, while 
radial orbits of infalling galaxies enhance the decline. 
The integral a{R) profile increases gradually to its max- 
imum value at ~ 1 Mpc, then steadily decreases before 
flattening out at a radius ~ 2.5 Mpc, indicating the a 
measured at this radius is no longer affected by velocity 
anisotropy, and is representative of the cluster potential 
(iFadda et al.l[l99l . 

We include the estimate of the velocity dispersion de- 
rived from the cluster X-ray temperature measured in 
^ ^2.3.11 bv assuming f3 = ^uipa'^/kT = 1 (i.e., assuming 
the same specific kinetic energy in the gas and galaxies), 
where a is the galaxy velocity dispersion, fi is the mean 
molecular weight of the gas particles, rrip is the proton 
mass and kT is the gas temperature within R ~ 500 kpc. 
The integral and differential ct{R) profiles appear, to 
within the errors, to be consistent with the value es- 
timated from the X-ray temperature when considering 
regions within 1 Mpc. 

3.3.3. 3-D tests for sub-structure 
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Fig. 12. — Differential velocity dispersion profile {top panel) 
and mean peculiar velocity (bottom panel) profiles as a function 
of projected radius. The solid lines show the total cluster values, 
while the dashed lines and error bars show the Icr confidence limits 
derived using the jackknife technique. The dot-dashed line shows 
the velocity dispersion derived from the mean X-ray temperature 
(see text), with the dotted lines showing the upper and lower limits 
based on the errors in the temperature (see § 12. 3. Il l 

The most efficient way to detect real physical substruc- 
tures is to use a combination of the 2-D spatial and 1-D 
velocity information t o search for local var iations in the 
velocity distribution (jPinknev et al.l [l996f ). The most 
common method h ere is to util i ze the A statistic de- 
veloped by Dressie r fc ShectmanI ()1988[ ). which tests for 
differences in the local mean and dispersion compared 
to the global mean and dispersion. The downfall of this 
method is that it assumes the global and local peculiar 
velocity distributions are Gaussian, which may not be 
true for dynamically acti ve syste ms. Here we pr efer the 
k-statistic, k, employed bv lCoUess fc Dunni (jl996f ). which 
is very similar to the A statistic but does not require the 
assumption of Gaussianity to be made. To determine k, 
the n = y/N nearest neighbors in projection are selected 
for each cluster member, where N is the total cluster 
member sample size, and the local velocity distribution 
of the n nearest neighbors is compared to the global clus- 
ter velocity distribution (minus the n nearest-neighbor 
velocities). Local departures from the global velocity dis- 
tribution are quantified using the K-S D statistic, with 
the null hypothesis that the local distribution is drawn 
from the global one and the significance determined by 
measuring the probability that the D statistic is larger 
than the observed D statistic for the observed sample 
size, Pks{.D > Dobs)- Then, k„ is defined as 

TV 

= ^ -\ogPKs{D > Dobs), (7) 
1=1 

giving a global measure of the substructure present in the 
cluster by summing the individual n values. The signifi- 
cance of Kn was determined by performing 10,000 Monte 
Carlo realizations with the peculiar velocities randomly 
shuffled, whilst maintaining the positional information, 
and remeasuring k„. 

The observed value of k„ = 379, is larger than any 
value obtained in our 10,000 realizations, which fol- 
low a log-normal distribution with a mean value of 
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Fig. 13. — Bubble plot outputs from the k test. The bold bub- 
bles are those deemed to be significant insomuch as they only occur 
in 5% of 10,000 realizations. The blue bubbles are those galaxies 
which have negative fpec, and the red bubbles have positive i;pec- 
The contours are galaxy density contours generated from applying 
a variable width Gaussian filter to the spatial distribution of spec- 
troscopically confirmed members. The contours are linearly spaced 
by 10 in the interval 10-150 galsMpc"^. Note the four dumpings 
of significant bubbles are coincident with overdensities in the pro- 
jected galaxy density. The center of the cluster is located at 0,0 
Mpc. 

^{InKsim) — 4.9 and standard deviation aihiKsim) = 
0.2. Thus, the observed k„ Hes 6(7 from the mean 
of the reahzations, and we can put an upper hmit on 
the probabihty of observing this k„ value by chance 
at less than 10~^. We conclude, therefore, that there 
is velocity substructure present in Abell 1201 at high 
significance. The results of the k test are best pre- 
sented using "bubble plots" showing a circle with radius 
r cx ~\ogPKs{D > Dobs) at each galaxy position, such 
that clustered large bubbles reveal local departures from 
the global velocity distribution. We show these bubbles 
in Figure [T3] and color code them based on the sign of 
the peculiar velocity, with blue and red having negative 
and positive Wpec, respectively. Overplotted are contours 
of galaxy surface density which have been produced by 
applying a variable width Gaussian filter, with a varying 
from 100 kpc in the cluster center to ^ 400 kpc in the 
outskirts, to the spatial distribution of the spectroscop- 
ically confirmed members. We define significant values 
of — log Pks{D > Dobs) as those which occur only 5% 
of the time in the 10,000 Monte Carlo realizations, and 
these are highlighted by the bold bubbles in Figure [T51 
Visual inspection of the distribution of large clustered 
bubbles in Figure [13] confirms the significance of the mea- 
sured Kn, as there appear to be four conglomerations of 
significant bubbles, each coincident with an increase in 
projected galaxy density. 

Having verified the existence of substructure within 
Abell 1201, we further utilize the spatial and velocity 
information of our cluster galaxies by allocating group 
membership using the Kaye's Mixture Model (KMM) al- 



gorithm of I Ashman et al.l ()1994[ ). The algorithm fits a 
user specified number of N-dimensional Gaussians to the 
data and determines the improvement of the fit over that 
of a single N-dimensional Gaussian via a maximum likeli- 
hood test. The major drawbacks of this method are that 
the spatial distribution of the galaxies does not follow a 
Gaussian shape (although the velocity distribution does 
for a relaxed cluster), and the number of Gaussians needs 
to be known a priori. Visual inspection of the projected 
X and Y galaxy distributions reveal they are at least 
qualitatively Gaussian, and given the benefit of includ- 
ing an extra two dimensions in the analysis far outweighs 
the false assumption of Gaussianity, we proceed with the 
full 3-D KMM analysis. Overcoming the latter drawback 
requires a robust method of estimating the initial number 
of Gaussians, and also their parameters. Given the corre- 
lation between the significant bubbles and galaxy surface 
density peaks seen in Figure[T3l we use it as a guide to es- 
timate the positions and projected radii of substructures. 
There appear to be 6 spatially separated substructures 
where the local velocity distribution differs significantly 
from the global one, along with the main Abell 1201 clus- 
ter. We inspect the velocity distributions of all galaxies 
within the estimated projected radius for each substruc- 
ture, excise any obvious interlopers and determine the 
median and standard deviation of the X position (kpc), 
Y position (kpc) and velocity distributions for the re- 
maining substructure galaxies. These parameters serve 
as initial estimates for input into the KMM algorithm, 
and are presented in tabled along with the outputs from 
the KMM algorithm, where {x, y, v) are the means of the 
distributions, {cx, (Jy, ay) the standard deviations, Ngai 
is the number of galaxies in the substructure and Rate is 
the estimate of the overall rate for correct allocation of 
galaxies to this substructure. 

Given that the substructures KMM3 and KMM4 are 
close both spatially and in velocity, as are KMM5 and 
KMM6, it is possible that they are part of the same struc- 
tures. We therefore combined the inputs for the four sub- 
structures into KMM(3-h4) and KMM(5-H6) and re-ran 
the KMM algorithm on the 5 substructures. The results 
for KMMl, KMM2 and KMM7 were very similar to those 
found using 7 partitions, with similar galaxies being allo- 
cated to the combined KMM(3-H4) and KMM(5-|-6) sys- 
tems as were allocated when considering them as sepa- 
rate entities, and the overall correct allocation estimator 
for both the 5 and 7 substructure analyses being 98% and 
97.9%, respectively. It appears the algorithm does not fa- 
vor a 7 structure partition over a 5 structure partition, 
so we used the simpler decomposition and assumed the 5 
substructure partition is the correct one. We tested the 
assumption of Gaussianity for the velocity distribution 
of each substructure using a K-S test and found no ev- 
idence for significant deviations from Gaussianity. The 
results of the allocations to the different substructures 
are plotted in the top panel of Figure [Ml where we have 
color-coded the symbols to match the respective velocity 
distributions plotted in the bottom panel of Figure [141 

4. MERGER SCENARIO 

Both the X-ray and optical analyses give clear indi- 
cations that Abell 1201 hosts multiple substructures. 
In this section we first give a qualitative scenario ex- 
plaining the appearance of the X-ray structure and cold 
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TABLE 4 

Results of the KMM analysis for partitioning the data into 7 (top) and 5 (bottom) substructures. See 
Figure [14] FOR the spatial and velocity distributions of the 5 partition results. 





Initial Inputs 






KMM Outputs 






Group 


(x, y,v) 






(x,y,v) 


(ca; , (Ty , (Ti, ) 




Rate (%) 


7 substructures 


KMMl 


(138, 388, 444) 


(103, 64, 175) 


9 


(150, 383, 432) 


(145, 121, 167) 


12 


100 


KMM2 


(-1125, 550, -1816) 


(320, 283, 318) 


11 


( -890, 346, -1887) 


(419, 224, 254) 


14 


100 


KMM3 


(-431,-444, -996) 


(62,134, 221) 


7 


( -462, -427, -1008) 


(44, 139, 160) 


6 


100 


KMM4 


(-869,-63, -607) 


(134, 116, 383) 


12 


( -924, -58, -672) 


(146, 120, 257) 


12 


100 


KMM5 


(1563, 513,-642) 


(200, 139, 233) 


11 


(1650, 533, -769) 


(180, 132, 208) 


10 


98 


KMM6 


(2150, 1213, -236) 


(155, 249, 170) 


9 


(2037, 1429, -317) 


(401, 459, 220) 


15 


96 


KMM7 


(0, 0, 2) 


(1236, 979, 738) 


262 


( -96, 93, 133) 


(1211, 1016, 682) 


252 


98 


5 substructures 


KMMl 


(138 388, 444) 


(103, 64, 175) 


9 


(151, 382, 432) 


(145, 121, 166) 


12 


100 


KMM2 


(-1125, 550, -1816) 


(320, 283, 318) 


11 


(-890, 346, -1886) 


(419, 224, 255) 


14 


100 


KMM(3+4) 


(-720, -200, -777) 


(239,223, 321) 


19 


(-684, -249, -867) 


(304, 264, 211)) 


16 


94 


KMM(5+6) 


(1884, 849, -608) 


(323, 421, 246) 


20 


(1831, 900, -753) 


(348, 452, 328) 


17 


99 


KMM7 


(0, 0, 2) 


(1236, 979, 738) 


262 


(-32, 138, 127) 


(1250, 1028, 663) 


262 


98 


Note. — The units of x, y, and 


(Ty arc kpc, and the units of 


V and (7^ arc kms^^ 









fronts through simple interpretations of the optical and 
X-ray o bservations combin ed wi th hydrodynamic simula- 
tions of [Poole et all (|2006f) and I Ascasibar fc MarkevitchI 
(|2006l ). Second we use two-body analytic models to de- 
termine which of the substructures are bound to the main 
cluster. 

4.1. KMMl, the north-west X-ray excess and the 
formation of the cold fronts 

Figure [2] shows an SDSS r-band image of the central re- 
gions of Abell 1201 with X-ray contours overlaid, along 
with regions showing the cluster members and KMMl 
allocations. Clearly, KMMl is coincident with the ex- 
cess X-ray emission which is probably the remnant of 
the gas core of KMMl. The morphology of the X-ray 
emission indicates the remnant KMMl core is breaking 
up, and there appears to be a tail pointing towards the 
main cluster, although it is difhcult to disentangle the 
cluster and subclump emission. The positioning of the 
cold fronts on opposite sides of the cluster center, along 
the direction to the north-west clump, suggests motion 
of the cluster core in this direction. This evidence points 
to a scenario where KMMl has made its closest approach 
to the core of Abell 1201 and is traveling outwards to- 
ward the north-west. Projection effects make it difficult 
to know on exactly which plane the merger is occurring, 
although the small radial velocity offset between KMMl 
and Abell 1201 and the location of the cold fronts sug- 
gest the majority of the subcluster motion is in the plane 
of the sky. Since the core of Abell 1201 is not completely 
disrupted and is coincident with the dominant cluster 
galaxy, which presumably lies at the cluster potential 
minimum, it is u nlikely KMMl pas sed directly through 
the cluster core (|Pooie et all l2008f ). The low velocity 
dispersion and compact galaxy distribution of KMMl 
suggest we are seeing the surviving central region of a 
once larger cluster which has been stripped of its outer 
members due to the tidal effects of the main Abell 1201 
cluster potential. 

The south-east cold front appears to be connected to 
the core of the main cluster, as evidenced by the signifi- 



cant residuals extending from the core to the cold front, 
and also by the finger of cold (~ 4 keV) gas joining the 
core and the cold front seen in the temperature map. 
This suggests the gas causing the cold front was once 
part of the core and has been displaced by the merger. 
It is possible that the cold fronts we see are the result 
of spirals of intertwined low and high entropy gas, as 
seen in both simu lations (jAscasibar &: Markevitchll2006 

; all 



Poole et al] I2006D and observations ([Fabian et 



2006 



Clarke et al.l 120041) . with the rotation axis of the spirals 



roughly in the plane of the sky. This scenario is capa- 
ble of producing several features observed in Abell 1201: 
two cold fronts at different distances on opposite sides 
of the core, the coincidence of the X-ray core with the 
dominant cluster galaxy and the region of hot gas on 
the north-west side of the cluster core which was prob- 
ably heated to its observed temperature by compression 
or shock heating caused by the merger with KMMl. To 
illustrate this we present Figure [12] which shows three 
temperature maps from the offse t (r^m — 360 kpc) , 3 to 
1 mass ratio merger simulation of lPoole et al.l ([20061 ) gen- 
erated using emission weighted temperatures integrated 
along a 3.5 Mpc line-of-sight. There are 3 projections 
each at 0.7Gyrs after pericentric passage with one pro- 
jection viewing along an axis perpendicular to the plane 
of the merger orbit {x — y) and two projections orthogo- 
nal to this. This makes apparent the effect of projection 
- the spiral type structure produces cold fronts observ- 
able over a wide range of viewing angles. Qualitatively, 
the map shown in the middle panel of Figure [15] repro- 
duces the features of Abell 1201 noted above, although 
the secondary core appears farther from the cluster cen- 
ter than that in Abell 1201. If the orbit of the secondary 
core is close to our line-of-sight and the core is close to 
turn around, this might be resolved by projecton effects. 
However, the simulation illustrated has generic initial 
conditions, not specific to Abell 1201. In particular, it 
is possible that a closer core passage could produce the 
structure observed in Abell 1201 without appealing to 
fortuitous viewing conditions. The fact that an ideal- 
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Fig. 14. — Top Panel: The spatial distribution of the different 
partitions assigned by the KMM analysis. The ellipses show the 
2(T contours for the Gaussians fitted to the respective spatial distri- 
butions. We also plot the galaxy surface density contours in pink 
(same spacing as Fig. Bottom panel: The velocity distribu- 

tions of the different partitions; the curved lines represent Gaussian 
functions whose mean and a are equal to the KMM v and Ct, val- 
ues. The color coding in the bottom panel matches the key in the 
top left of the top panel. The bottom right-most velocity distri- 
bution shows the whole cluster sample with a Gaussian generated 
using the biweight estimators overplotted in pink; the combination 
of all the KMM Gaussians is plotted in black. 

ized simulation reproduces the majority of the observed 
features supports the proposed merger scenario. Future 
simulations specific to Abell 1201 will determine which 
of these possibilities can better explain its observed fea- 
tures. 

It appears Abell 1201 is an excellent example of a 
"sloshing" type cold front cluster, where the perturber is 
still clearly visible in the X-ray and optical observations 
(i n the form of KMMl), similar t o Abell 1644 (see Figure 
IT lMarkevitch fc Vi khlinin"2007[). Thus, Abell 1201 is an 
excellent candidate for follow up detailed simulations in 
order to derive an accurate picture of exactly how the sys- 
tem has evolved to its current state, and whether it will 
evolve into a relaxed looking cluster harboring cold fronts 
with no discernable perturber (eg. , RXJ1720. 1-1-2638, 
MS1455. 0+2232 or Abell 2029 Markevitch fc VikhlininI 
I2007t iMazzotta fc Giacintuccill2008f i 



4.2. Two-body merger dynamics 

In this section we appl y the two-body dyn amical anal- 
ysis first implemented by 'B eers et al.l ()1982[ ) to the sub- 
structures detected in §3. 3. 31 with the aim of determining 
which of the substructures are bound to the main clus- 
ter, and to give an initial idea of the internal dynamics 
of the cluster. The model allows estimation of the prob- 
ability that a substructure is unbound and lying close to 
the main cluster along the line of sight, and also allows 
the calculation of the probability that a substructure is 
bound and collapsing or bound and expanding. 

The accuracy of the results given by the model rely 
critically on the following assumptions: the orbits are 
radial, the cluster masses are concentrated into a point 
at the respective centers, the clusters had zero initial 
separation at t=0, and they are moving apart or com- 
ing together for the first time. These assumptions are 
extremely unlikely to hold true for KMMl, which prob- 
ably resides well within the cluster virial radius, deep 
within the main cluster potential well (where dynamical 
friction, tidal forces and angular momentum all become 
significant) and is apparently currently heading to the 
north-west after pericentric passage. Thus, we do not 
present results from the two-body analysis for KMMl, 
and defer quantitative analysis of the KMMl merger un- 
til detailed simulations are available, noting only that 
given the observations presented here, it appears KMMl 
and Abell 1201 form a bound system. 

As inputs, the model requires information about the 
projected spatial separation, Rp, the line of sight velocity 
difference, Vr, and the total mass of the system, returning 
possible solutions for a (the angle between the line join- 
ing the two clusters and the line of sight), the total mass 
required to bind the system, and the true 3-D spatial 
separation, R, and velocity difference, V. The paramet- 
ric solutions to the equations of motion for bound radial 
orbits are 



V = 



sina 



/2GA/y/^ sinx 
V Rm J (1 - cosx) ' 



^ ( \ 

\8GM J 



1/2 



(X - sinx) , 



R = = ^ (1 - cosx) , 
cosa 2 



(8) 



(9) 



(10) 



where Rm is the maximum separation of the two clusters 
at turn around, M is the total mass of the system, and x 
is the developmental angle which varies from < x < 27r 
with X = 0, 27r being the stages of the orbit corresponding 
to zero spatial separation. It is also possible to solve for 
the unbound case, where the parametric solutions are 



V 



Vr 



sma 



sinhx 



(coshx — 1) ' 



R = 



GM , . , 

t = TTT (smhx - X) , 

Rp GM , , 
- -7T (coshx - 1) 



cosa 

where Voo is the asymptotic expansion velocity. 



(11) 
(12) 

(13) 
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Fig. 15. — Snapshot temperature maps of the 3 to 1 mass ratio, offset {r^i„ = 360 kpc) merger simulation from'Poolo et al.' f2006') taken 
0.7Gyrs after pericentric passage. The left panel shows the merger viewed from a vantage point where the line of sight is perpendicular 
to the merger orbital plane, while for the middle panel the line of sight is along the x-axis and the right panel the line of sight is along 
the y-axis. The vectors on the lower left of each panel show this, and the length of each vector has physical size 500 kpc. The color scale 
runs from black to white with black showing the lowest temeperatures, and white the highest (on an arbitrary temperature scale). X-ray 
surface brightness contours are overlaid and the white crosses show the positions of the peak in the dark matter density for the primary 
and secondary clusters. 

For the bound case, combining equations [51 [D and [TUl 
gives 

V,.t (1 - cosx)^ 



derived for galaxies assigned to KMM7 and 



tana 



Rp sinxix - sinx) ' 



and, similarly, we obtain 



tana 



Vi-t (coshx - 1)^ 
Rp sinhx(sinhx - x) 



(14) 



(15) 



Rpv 



■^viri^-^vir 1) 



for the unbound case. Assuming t — 11.38 Gyrs, i.e. the 
age of the Universe at z = 0.168 in our assumed cosmol- 
ogy, for values of the parameter x in the range to 2n 
these equations determine possible solutions for a. Given 
that the mass of the system is the least well constrained 
of the input parameters, we solved for M as a function of 
a using as input Vr and Rp for each group in the KMM 
analysis. Here Rp was calculated as the distance from 
the KMM centroid to the central dominant galaxy, rather 
than the KMM7 center, which is slightly offset from the 
central dominant galaxy position. Vr is the radial veloc- 
ity offset of the KMM partition of interest, taken with 
respect to the biweight location of peculiar velocities in 
KMM7, where either the median (for KMM groups with 
Ngai < 15) or the biweight locator {Ngai > 15) was used 
to determine the group's peculiar velocity. We plot the 
possible solutions of M as a function of a in Figure [TBI 
Determination of solutions for the orbits of the systems 
requires knowledge of the total mass of the system. Given 
the velocity dispersions of the substructures in Abell 1201 
are much smaller than that of the main cluster, it is 
reasonable to assume that when compared to the mass 
of the main Abell 1201 cluster the contribution of the 
masses of the subclumps is negligible. Hence we assumed 
the total system mass was equivalent to t he virial mass 
of Ab ell 1201. We used the methodology of lGirardi et al.l 
(|1998l) for the determination of the mass 



Stt afjRpv 
~2 G~ 



C, 



(16) 



where C is the surface term correction accounting for the 
lack of coverage for the entire cluster (i.e. out to the turn 
around radius) . Here we took C = 0.19M„,>, which is the 
median value derived in [Girardi et al.[ (|1998[ ). and cr^ = 
665 ± 32kms~^ which is the biweight scale estimator 



1 p-l' 



(17) 



which w as determined within r2oo — OAla^/ H{z) = 
1.5 Mpc ([Carlberg et al.[[1997f) where = 154 is the 
number of KMM7 members within r2oo, and Rij is the 
projected distance between two galaxies. We determined 
Rpv = 1.5 ± 0.1 Mpc and M = 5.9(±0.6) x lO^-^Mo 
where the error in Rpv was determined using the Jack- 
knife technique, and standard error propagation was used 
to derive the error in the mass. 

The possible solutions for the orbits are given by the 
intersection of the line representing the Abell 1201 virial 
mass estimate with the curve for M as a function of a in 
Figure [inj For each solution, we present the associated 
probability relative to the other solutions by considering 
the range of possible a given by the upper and lower 
values, au and a^, respectively, derived from the error 
bounds on each curve, and assuming that each individual 
solution is equally probable. The relative probability is 



Prel 



( cosada 



for the solution of interest, and 



Prel = PBO + PUO + PBI^ + PBh 



(18) 



(19) 



where pbo, Puo, Pbi^, Pbi^ are the J^^" cosada values 
for the bound outgoing, unbound outgoing, first bound 
incomi ng and second boun d incoming solutions, respec- 
tively ([Brough et al.l I2OO60 . We present the results in 
Table [Sj 

We also determined the probability the subclumps are 
bound using the Newtonian binding criterion, where a 
two-body system is bound if the kinetic energy is less 
than or equal to the poten tial energy. The criterion is 
given by [Beers et al.l ([1982( ) as 



V^Rr, < 2GM sin^acosa. 



(20) 



This curve is also plotted in Figure [THl and the prob- 
ability is Pbound = J^" cosada, where a^ and au are 
determined from the intersections of Myir and the curve 
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Fig. 16. — Binding mass as a function of a {orange and green curves). The orange curves correspond to bound solutions, while the green 
curves represent unbound solutions. The dashed curves show the Icr errors. The blue line shows the mass derived using galaxies allocated 
to KMM7 and the dashed blue lines indicate the upper and lower Icr errors. The small black dashed line delineates the bound and unbound 
regions according to the Newtonian binding criterion. 



TABLE 5 

Probabilities that the detected substructures are bound to the Abell 1201 cluster, 
derived using the two-body dynamical analysis and the newtonian criterion (see text). 



Subclump 


Vr 


Rp 


Pbi^ 




Pbo 


Puo 


Pbound 




{ km s^-^) 


(Mpc) 


(per cent) 


(per cent) 


(per cent) 


(per cent) 


(per cent) 


KMM2 


-1989 ±46 


0.95 








100 




KMM(3+4) 


-996 ± 72 


0.73 


62 


38 




« 1 


60 


KMM(5+6) 


-798 ± 113 


2.04 






« 1 


100 


41 



representing the Newtonian criterion. We present the 
results in Tabled 

Considering the Newtonian criterion and two-body 
sohitions, we find that there is a high probability of 
KMM(3-h4) being bound to AbeU 1201, with two solu- 
tions which are, roughly speaking, as probable as each 
other. The most probable solution is KMM(3-|-4) is 
bound and incoming with a = 26 degrees, R — 0.8 Mpc 
and V = —2252 km s^^, whilst the slightly less probable 
solution is also bound and incoming with a = 71 degrees, 



R = 2.3 Mpc and V = -1052 km s^^. If KMMl is out- 
bound to the north-west, as proposed in § 14.11 then 
it is possible that KMM(3+4) consists of tidal debris 
stripped from KMMl during is approach to pericentric 
passage from the south-east. Further evidence for a tidal 
origin comes from the spatial distribution of galaxies, 
where KMM(3-f4) appears as two conglomerations in 
Figure [TBI although simulations are required to confirm 
this. 

In the case of KMM2, we find only unbound out- 
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going solutions with the most probable solution giving 
a — —88 degrees (i.e. it's motion is aligned at 2 deg to 
the line of sight), R = 24.8 Mpc and V = 1990 km s'^ 
Thus it is likely that KMM2 is a group of galaxies ly- 
ing in the foreground moving away from Abell 1201, and 
not physically associated. KMM(5-|-6) has a best fit- 
ting unbound outgoing solution, although within the er- 
rors this may be bound outgoing, with a — —81 degrees, 
R = 12.9 Mpc and V = 808 km s^^. Given the Newto- 
nian criterion allows KMM(5-f 6) to be bound and that it 
is possible we have underestimated the mass within the 
radius at which KMM(5-|-6) lies, it is entirely plausible 
that we have also underestimated the binding probability 
for KMM(5-h6). 

The two body analysis presented here has allowed the 
determination of which substructures are most likely 
bound to the main cluster. This gives an initial under- 
standing of which substructures are important in under- 
standing the internal dynamics of Abell 1201, which will 
be further modeled in future simulations and presented 
in a forthcoming paper. 

5. SUMMARY AND CONCLUSIONS 

We have presented an analysis of Abell 1201 using both 
Chandra X-ray data and optical spectroscopic data ob- 
tained using the AAT/AAOmega and MMT/Hectospec 
MOS. The X-ray analysis reveals the following: 

1. An elliptical morphology with two surface bright- 
ness discontinuities positioned roughly on an axis 
joining the main cluster with a clump of excess 
emission to the north-west. The density, temper- 
ature and pressure jumps across these discontinu- 
ities are consistent with cold fronts. 

2. A residual map obtained by subtracting a double 
beta model from the surface brightness distribution 
shows significant residuals at the position of the 
north-west clump. It shows significant residuals 
extending from the cluster core to the south-east 
cold front. 

3. A temperature map shows that south-east residual 
is coincident with a finger of cold gas which also 
extends from the cluster core and terminates on 
the inner side of the front. 

4. There is also a significantly hotter region on the 
north-west side of the inner cold front between the 
main cluster core and the subclump further to the 
north-west. 

The optical MOS analysis reveal: 

1. From 321 cluster member spectra the cluster red- 
shift is z=0.1673 and the velocity dispersion is 
778 km s-^ 

2. The velocity distribution is not significantly non- 
Gaussian, despite the clearly disturbed nature of 
Abell 1201 evidenced by the X-ray data. 

3. Combining the peculiar velocity and spatial infor- 
mation reveals significant localized velocity sub- 
structure. 



4. Using the KMM method of I Ashman et all ()1994f ) 
the substructure can be partitioned into 5 dis- 
tinct clumps - the main Abell 1201 cluster with 
262 members (KMM7), an infalling subgroup ^ 
410 kpc to the north-west coincident with the X- 
ray excess containing 12 members (KMMl), a sub- 
group ^ 950 kpc to the NE with 14 members 
(KMM2), a subgroup ~ 730 kpc to the south-east 
with 16 members (KMM3-f 4) and a group ~ 2 Mpc 
to the north-west with 15 members (KMM5-I-6). 

5. Application of the two-body dynamical analysis of 
iBeers et all (|l98a reveals KMM2 is likely an un- 
bound foreground group, whilst it is also possible 
KMM(5-h6) is unbound and KMM(3-h4) has a high 
probability of being bound to the main cluster. 

The analysis presented here supports the view that the 
cold fronts in Abell 1201 are a direct result of merger 
activity. Our X-r ay and optica l anal ysis combined with 
the simulation of iPoole et al.l ()2006f ) indicate that the 
KMMl subcluster has come from the south-east and is 
currently heading towards the north-west after an offset 
passage past the main cluster core. Disturbance caused 
by the passage of this subcluster has set off sloshing of the 
cool gas in the core of the main cluster, which produced 
two concentric cold fronts around the main core. The 
value of the use of multiple methods of merger detection 
for analyzing a clusters dynamical state is also confirmed 
in this study. 

The next step in our study of cold front clusters will 
be to utilize radio data in combination with the optical 
spectroscopy presented here to study the effects of the 
merger on the galaxy population, and search for signs 
of non-thermal radio halo/relic emission. Deeper X-ray 
observations will help place tighter constraints on the 
parameters derived here, whilst allowing a more detailed 
analysis of the complex X-ray structure. 
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